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ABSTRACT 


A  \ 

A  method  is  developed  to  analyze  and  compare  the  effectiveness  of 

ASW  vehicles*  The  measure  of  effectiveness  is  the  probability  that  the 
vehicle,  after  detecting  a  submarine  with  passive  sensors,  can  transit 
to  the  contact  area  and  re-establish  contact  with  the  submarine*  A  com¬ 
puter  simulation  is  developed  and  an  example  using  three  hypothetical 
ASW  vehicles  is  illustrated. 
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Introduction. 


As  a  result  of  increased  emphasis  currently  being  placed  on  ASW 
weapon  systems,  it  is  advantageous  to  possess  a  computer  simulation 
capable  of  analyzing  and  comparing  proposed  ASW  vehicles  and  opera¬ 
tional  ASW  vehicles.  The  computer  simulation  for  the  comparison  of  ASW 
vehicles  developed  in  this  paper  is  intended  to  be  used  for  this  pur¬ 
pose. 

The  problem  at  hand  is  to  develop  a  method  for  comparing  ASW  vehi¬ 
cles  under  operational  conditions  using  vehicle,  submarine,  and  tactical 
situation  parameters*  The  simulation  allows  one  to  compute  a  proba¬ 
bility  which  is  used  as  a  measure  of  effectiveness*  This  probability, 
hereafter  referred  to  as  the  "probability  of  success,"  is  defined  in 
this  paper  as  the  probability  that  an  ASW  vehicle  successfully  relocates 
a  submarine  at  an  estimated  position  given  that  a  detection  has  occurred, 

A  description  of  the  simulation  and  its  Imitations  is  made  in  Sec¬ 
tion  2*  It  is  very  important  to  understand  the  simulation  and  the 
background  assumptions  before  the  results  of  the  simulation  are  ana¬ 
lyzed*  After  the  problem  is  developed  mathematically  in  Section  2,  a 
discussion  and  an  illustration  is  given  in  Section  3  to  demonstrate  how 
to  use  the  computer  program.  In  Section  4,  several  applications  are 
discussed*  Finally  a  comparison  is  made  of  three  hypothetical  ASW 
vehicles- -a  High  Speed  Surface  Vehicle,  a  VERTOL  Vehicle,  and  a  Seaplane- 
Type  Vehicle- -using  the  results  generated  from  the  computer  simulation. 
The  output  of  one  computer  run  is  illustrated  in  Figures  5  and  6.  Curves 
for  the  probability  of  success  versus  range  to  the  submarine  for  each 
vehicle  are  displayed  in  Figures  7,  8,  and  9*  These  curves  are  used  to 
analyze  and  compare  the  systems* 


1 


The  results  of  these  comparisons  Illustrate  the  usefulness  of  the 
simulation  in  comparing  the  performance  of  different  vehicles  and  for 
studying  the  effect  of  various  parameters  such  as  expected  sensor  range 
on  the  probability  of  success* 
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2.  Simulation  Description  and  Mathematical  Development. 

Effort  has  been  directed  toward  the  development  of  a  computer  simu¬ 
lation  that  represents  a  real  world  conflict  between  an  ASW  vehicle  and 
a  submarine.  The  setting  is  one  of  an  ASW  vehicle  assigned  to  a  patrol, 
surveillance,  or  barrier  mission.  Initial  detection  is  assumed  to  have 
been  made  using  passive  sensors.  The  vehicle  then  proceeds  as  rapidly 
as  possible  to  an  estimated  position  (EP)  in  order  to  relocate  the 
submarine.  The  EP  is  the  predicted  position  of  the  submarine  when  the 
actual  position  is  not  known,  but  the  vehicle* s  sensors  indicate  the 
presence  of  a  submarine.  The  objective  of  the  vehicle  is  to  compute  an 
EP  of  the  submarine,  close  the  position  at  maximum  speed,  and  attempt 
to  relocate  the  submarine  using  active  sensors.  No  information  about 
the  actions  of  the  submarine  is  available  while  the  vehicle  closes  to 
the  EP. 

Certain  other  restrictions  are  assumed  in  the  model.  Ahead  thrown 
weapons  and  nuclear  weapons  are  not  available  for  use  by  the  vehicle. 
Multiple  targets  are  excluded.  A  passive  detection  range  of  at  least 
five  nautical  miles  is  necessary;  otherwise,  an  active  sensor  would  be 
used  and  the  problem  of  transiting  to  an  estimated  position  is  irrele¬ 
vant.  Also,  the  vehicle  is  assumed  to  lie  motionless  in  the  water  until 
detection  has  been  achieved  to  insure  passive  sensor  capability. 

The  simulation  is  developed  for  two  different  tactical  situations. 
The  first  case  is  that  of  a  submarine  which  is  unaware  of  the  activity 
of  the  vehicle.  The  second  case  is  that  of  a  submarine  which  is  alerted 
by  the  activity  of  the  vehicle.  Common  to  both  of  these  cases  is  a 
maneuvering  board  problem  which  will  be  discussed  first. 
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The  maneuvering  board  solution  consists  of  describing  geotoe  trie  ally 
the  positions  and  motions  of  a  vehicle  and  a  submarine  with  time.  The 
passive  sensor  of  the  vehicle  provides  the  last  known  position*  course 
and  speed  of  the  submarine*  and  the  bearing  and  range  to  the  submarine 
from  the  vehicle.  Knowing  these  values  and  the  vehicle  speed*  an  EP, 
transit  time  of  the  vehicle,  relative  closing  speed,  and  the  distance 
traveled  to  EP  are  computed*  For  overall  ease  in  reading  the  same  nota¬ 
tion  is  used  in  the  mathematical  development  that  has  been  used  in  the 
computer  simulation  (Appendix  C). 

The  solution  is  computed  by  solving  a  velocity  vector  problem*  The 
vehicle  is  located  at  the  origin  (0,0)  of  a  rectangular  coordinate  sys¬ 
tem,  One  side  of  the  velocity  triangle  (Figure  1)  represents  the 
estimated  velocity  vector  of  the  submarine  (Ul)  which  is  re-established 
at  the  origin.  The  velocity  vector  for  the  vehicle  (E)  is  unknown  in 
direction,  but  known  in  magnitude.  The  remaining  relative  velocity  vec¬ 
tor  (RELSPD)  represents  the  relative  course  and  speed  of  the  vehicle 
closing  the  submarine.  Only  the  direction  of  RELSPD  is  known  which  is 
parallel  to  the  true  bearing  line  from  the  vehicle  at  (0,0)  to  the  sub¬ 
marine  at  (X2,  Y2).  Then  the  intersection  (X3 ,  Y3)  between  the  relative 
velocity  vector  and  the  vehicle  velocity  circle  is  computed. 

In  order  to  solve  the  maneuvering  board  problem  various  points  and 
vectors  in  the  rectangular  coordinate  system  must  be  computed*  The 
values  for  the  components  of  Ul  are; 


XI  =  Ul  x  SIN 

(THETA) 

(1) 

Y1  -  Ul  x  COS 

(THETA) 

(2) 

The  last  known  position  (X2,y2)  of  the  submarine  is  given  by: 
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X2  =  RAI  x  SIN  (BETA) 

(3) 

Y2  =  RAI  x  COS  (BETA) 

(A) 

where  BETA  is  the  bearing  in  degrees  of  the  submarine  from  the  vehicle 
and  RAI  is  the  range  in  nautical  miles  of  the  submarine  from  the  vehicle. 
The  RELSPD  vector  is  represented  by  the  line  through  the  point  (XI,  Yl) 
and  parallel  to  the  line  connecting  the  origin  and  (X2,  Y2).  The  equa¬ 
tion  of  this  line  is: 

Y3  =  ||  (X3-X1)  +  Yl  .  (5) 

The  RELSPD  vector  intersects  the  vehicle  velocity  circle  at  the  point 
(X3,  Y3).  The  equation  of  the  vehicle  velocity  circle  is; 

Y32  ■+■  X32  *  V2?  (6) 

where  V  is  the  average  vehicle  transit  speed.  The  point  (X3,  Y3)  is  de¬ 
termined  by  solving  equations  (5)  and  (6)  simultaneously.  The  direction 
of  the  vector  from  the  origin  to  (X3S  Y3)  represents  the  course  of  the 
vehicle  when  traveling  inbound  to  EP.  The  relative  speed  vector  is  rep¬ 
resented  by  the  vector  between  (XL,  Yl)  and  (X3,  Y3)  and  its  magnitude  is 

RELSPD  =  -\j  (X3-Xl)2+  (Y3-Y1)2.  (7) 

The  vehicle’s  transit  time  (TD)  to  EP  is  defined  as: 

TD  =  RAl/RELSPD  #  (8) 

Let  the  vehicle  reaction  time  (RT)  be  the!  following  sum: 

RT  =  PTOT  4-  TOT  +  CLOT  +  RELT  (9) 

where  y  PTOT  -  pre take-off  time  which  commences  when  the  vehicle  can 
no  longer  observe  the  submarine  and  includes  decision 
time,  warm-up  time,  and  sensor  retrieval  time, 

TOT  -  take-off  time  which  includes  take-off  roll  and  time  to 
take-off  into  the  wind. 


6 


CLOT  -  land  time  which  includes  time  to  land  and  taxi  time 

needed  to  arrive  at  EP  in  order  to  deploy  active  sensors, 

RELT  ~  target  relocation  time  which  includes  the  time  necessary 
to  deploy  sensors  and  ends  when  the  submarine  is  re¬ 
located. 

Any  of  these  times  can  have  a  zero  value.  For  example,  a  helicopter 
could  have  TOT  equal  to  zero  if  it  needs  no  take-off  time.  Also,  the 
times  are  mutually  exclusive  so  that  no  time  interval  is  added  twice. 
Obviously  the  blind  time  (TB)  is: 

TR  *  RT  4-  TD  ,  (10) 

Evasion  for  the  alerted  submarine  begins  when  blind  time  commences.  At 
that  time  the  submarine  is  assumed  to  be  at  (X2,  Y2).  If  RT  is  equal  to 
zero,  the  submarine  is  at  (X2t  Y2)  when  evasion  begins;  but  if  RT  is 
not  equal  to  zero,  the  submarine  is  not  actually  at  (X2,  Y2)  when  eva¬ 
sion  begins  because  it  takes  RT  for  the  vehicle  to  commence  moving 
towards  EP-  During  the  reaction  time  the  vehicle  Is  motionless,  but 
the  submarine  is  moving  away  from  (X2,  Y2).  Because  of  the  small  dis¬ 
placement  of  the  submarine  from  (X2,  Y2)  due  to  RT,  the  assumption  that 
the  submarine  is  at  (X2,  Y2)  when  evasion  begins  does  not  significantly 
affect  the  results  of  the  simulation. 

It  is  necessary  to  determine  the  estimated  position  (EP)  of  the 
submarine  after  a  time  interval  of  TB  and  a  submarine  speed  U3.  The 
equation  of  the  estimated  distance  (DIST)  traveled  by  the  submarine  is: 

DIST  =  U3  x  TB  .  (11) 

The  EP  is  also  the  position  at  which  the  vehicle  estimates  that  an  in¬ 
tercept  will  be  made  with  the  submarine.  The  components  of  EP  relative 
to  the  origin  are: 

X4  *  DIST  x  SIN  (THETA)  +  X2  (12) 
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Y4  *  DIST  x  COS  (THETA)  H-  Y2  * 


<U) 


The  equation  for  the  distance  DISDAT)  traveled  by  the  vehicle  to  EP  is 


(14) 


DISDAT 


The  values  for  TB  and  DIST  are  very  important  in  determining  the  proba¬ 
bilities  of  success  for  the  two  cases  developed  below* 

Case  I,  The  unalerted  submarine; 

It  is  conceivable  that  the  vehicle  could  transit  to  EP  without  the 
submarine  being  cognizant  of  the  impending  danger*  In  this  case  the  un- 
alerted  submarine  does  not  evade  and  proceeds  at  cruise  course  and  veloc¬ 
ity*  Prediction  errors  for  the  submarine  course  and  velocity  are  made 
which  are  characteristic  of  the  vehicle’s  equipment*  The  errors  are 
assumed  to  be  distributed  normally  about  the  predicted  course  and  veloc¬ 
ity,  Because  of  long  blind  times  expected,  the  predicted  positions  of 
the  submarine  are  not  accurately  represented  by  a  circular-normal  distri¬ 
bution,  and  a  numerical  integration  must  be  performed.  Probabilities  of 
success  are  computed  for  each  of  a  representative  distribution  of  pre¬ 
dicted  positions  consisting  of  25  points  or  cells  (see  Figure  2),  This 
concept  has  been  adopted  from  Q) j . 

The  last  known  submarine  position  (X2,  Y2)  is  established  at  the 
origin  of  a  rectangular  coordinate  system.  Since  DIST  is  the  distance 
traveled  by  the  submarine  during  TB,  then  EP  is  (DIST,D)  in  this  coor¬ 
dinate  system* 

Using  the  prediction  errors,  a  submarine  position  is  computed  for 
each  of  the  25  cell  midpoints*  The  distance  (VV)  and  angle  (THETER)  for 
each  cell  are  computed  using  the  submarine  cruise  speed  (Ul),  the  submarine 
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Figure*  2-  I^Vtgral'Io'A  FW  TW  Uttftlertexi  C^s4* 


i 


9 


speed  error  (UERR),  the  submarine  course  error  (THDEL),  and  the  blind 
time  (TB).  The  equations  for  W  and  THETER  are; 


W  =  (Ul-AA(I)  x  UERR)  TB  (15) 

THETER  =  AA  (J)  x  THDEL  ,  (16) 

with  AA(I)  and  AA(J)  s  0,  +  .52,  ±  1.28  depending  on  which  cell  is  being 
considered.  These  numbers  can  be  obtained  from  any  normal  density  tables, 
[2]. 

The  equations  for  the  submarine  position  (SUBX  (I,J),  SUBY(I,J))  at 
the  center  of  each  of  the  25  cells  are; 

SUBX  (I,J)  -  VV  x  COS  (THETER)  (17) 

SUBY  (I,J)  =  W  x  SIN  (THETER)  .  (18) 

A  bias  is  computed  between  EP  and  the  center  of  each  cell.  The  equation 
for  the  bias  of  the  (i,j)t^  cell  is: 

SBIAS  (I,J)  =  -\J  (SUBX(I,J)-DIST)2 (SUBY(I,J))2  .  (19) 

Various  errors  have  to  be  considered  to  give  the  simulation  realism. 
It  is  assumed  that  the  navigation  error  and  sensor  locating  error  are 
additive  about  the  point  EP.  The  distribution  of  the  vehicle  location 
about  EP  due  to  these  errors  is  assumed  to  be  circular  normal. 

The  standard  deviation  of  navigation  error  (SIGNAV)  is; 

SIGNAV  s  PRCNTN  x  RAI  (20) 

where  PRCNTN  is  a  percentage  of  the  last  known  range  (RAI).  The  stand¬ 
ard  deviation  of  sensor  locating  error  consists  of  the  standard  devia¬ 
tion  of  bearing  (  SIGSBR  )  and  the  standard  deviation  of  range  (SIGSR) 
which  are  defined  by: 

SIGSR  -  PRCNTS  x  RAI  (21) 

(22) 


SIGSBR  *  RAI  x  SIN  (BRGER) 
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where  BRGER  is  the  bearing  error  and  FRONTS  is  the  percentage  error  of 
the  last  known  range*  Since  (21)  is  greater  than  (22) ,  because  of  the 
passive  sensor,  an  elliptical  normal  distribution  exists;  but  it  is 
approximated  by  a  circular  normal  distribution  for  ease  of  computation. 
This  is  accomplished  by  equating  the  area  of  an  ellipse  to  the  area  of  a 
circle.  The  result  is  an  approximate  standard  deviation  (SXGAPR). 


7 rcr'  cr  slVo-'2  *  with 


XE 


YE 


CIR 


^~XE  =  SXGSBR 


0"ye  = 


SIGSR 


(23) 


SXGAPR 


0”cir  =~\J  CTxe  CTye 


(24) 


Therefore,  the  total  standard  deviation  of  error  (SIGJT)  is  the  sum 
of  the  squares  of  the  approximate  standard  deviation  of  error  and  the 
standard  deviation  of  navigation  error  according  to  the  equation: 


SIGJT  •~]j  (SIGNAV)2-*-  (SIGAPR)2  .  (25) 

Other  error  models  were  considered  in  order  to  more  accurately  represent 
the  errors  involved.  Due  to  lack  of  time  the  error  model  described 
above  *  although  not  the  best,  was  adopted. 

To  compute  the  probability  of  success  for  the  unalerted  submarine 
the  circular  coverage  function  is  usedf£lj,  In  order  to  use  this  func¬ 
tion,  the  parameters  SBIAS(I,J),  SIGJT,  and  EXPRNG  are  required.  The 
first  two  have  been  defined  above.  The  term  EXPRNG,  which  is  defined  as 
the  expected  range  of  the  active  sensor  used  by  the  vehicle,  is  a  charac¬ 
teristic  of  the  sensor  employed.  The  quantities  SB IAS (I , J)  and 

SIGJT 


EXPRNG 

SIGJT 


are  used  In  the  circular  coverage  function  to  compute  the  proba¬ 


bility  of  success  (PK(X,J))  that  a  circular  disk  of  radius  EXPRNG  will 
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cover  a  point  (SUBX(X,J),  SUBY(I,J)  from  the  EP  since  the  probable  posi¬ 
tion  of  the  disk  is  described  by  a  Gaussian  distribution.  This  proba¬ 
bility  PK(I*J)  is  the  probability  of  success  for  the  ijth  cell.  Since 
it  is  equally  likely  that  the  submarine  is  in  each  of  the  25  cells,  ♦,  "■ 
the  probability  that  it  is  within  each  cell  is  .04,  The  probabilities 
of  success  for  each  cell  are  summed  and  multiplied  by  .04  to  give  the 
"probability  of  success  (PROB)  against  an  unalerted  subtnar ine* 11  The 
probability  of  success  is; 


(26) 


PROB 


Case  II.  The  alerted  submarine: 

The  submarine  is  aware  of  the  presence  of  an  ASW  vehicle  and  evades 
by  changing  course  and  speed.  The  speed  changes  are  characteristics  of 
the  submarine  but  the  course  changes  are  determined  in  number  and  degree 
by  the  situation  being  simulated.  A  probability  of  success  is  computed 
for  each  evasion  turn  and  an  arithmetic  average  taken  over  the  number 
of  evasion  tactics  used. 

This  case  is  designed  to  give  practically  full  evasion  to  the  sub¬ 
marine  (see  Figure  3).  The  submarine  evasion  turns  are  limited  to  a 
maximum  turn  of  90  degrees  to  either  side  of  the  estimated  submarine 
course.  No  additional  turns  are  granted  after  the  initial  turn  is  exe¬ 
cuted,  For  turns  greater  than  90  degrees  the  submarine  would  deviate 
from  its  initial  track  to  such  an  extent  that  it  would  not  be  able  to 
make  up  the  lost  time  necessary  to  make  good  a  mission  speed  of  advance. 

The  range  of  values  for  the  set  of  evasion  turns  is  determined  by 
the  situation  being  simulated.  For  example,  a  slow  submarine  would 
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almost  always  turn  away  from  the  vehicle  when  the  vehicle's  presence  is 

known.  There  are  two  limiting  turns  which  bound  the  set  of  evasion  turns 

to  the  right  and  the  left.  The  bounds  must  be  no  greater  than  90  degrees 

to  either  side  of  the  estimated  submarine  course,  PSI  and  DELTA  are  the 

left  bound  and  right  bound,  respectively.  The  range  of  values  for  the 

set  of  evasion  turns  is  determined  by  (PSI-BELTA).  The  interval  (ETA) 

between  evasion  turns  is  ETA  s  PSI -DELTA  (  HU  being  the  number  of  eva- 

NU 

sion  turns  considered.  If  the  set  of  turns  encompasses  turns  to  the 

right  of  the  submarine  course  then  ETA  =  — (360-DELTA)  ^  |n  pigure  3 

NU 

PSI  =  060  degrees,  DELTA  =  330°,  and  ETA  =  060  (360-330)  r  15°. 

6 

Therefore,  the  evasion  turns  considered  for  this  example  are  330°,  345°, 
360° ,  015°,  030°,  045°,  and  060°. 

In  the  alerted  submarine  case  the  submarine  course  and  velocity  pre¬ 
diction  errors  that  were  present  in  the  unalerted  case  are  ignored. 

Since  in  this  case  a  set  of  evasion  turns  are  selected  which  are  only  sub¬ 
jective  predictions,  it  is  valid  to  exclude  prediction  errors  because 
they  are  comparatively  insignificant  to  the  computations.  However,  ve¬ 
hicle  sensor  and  navigation  errors  that  were  assumed  in  the  unalerted 
case  still  prevail. 

Initially,  the  submarine  is  at  cruise  speed.  When  evasion  confluences, 
the  assumed  submarine  speed  is  increased  to  fast  speed.  The  character¬ 
istics  of  the  submarine  determine  how  long  the  submarine  remains  at  fast 
speed  (length  of  time  is  Tl)  after  which  it  slows  to  silent  speed.  The 
time  intervals  of  acceleration  and  deceleration  are  considered  negligible 
when  compared  with  the  long  blind  times  encountered. 

The  alerted  case  uses  a  rectangular  coordinate  system  with  the  last 
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known  position  of  the  submarine  (X2,Y2)  at  the  point  of  origin  when  eva¬ 
sion  begins*  From  this  point  positions  (X(N) ,  Y(N))?  which  are  the 
positions  at  the  end  of  each  of  the  N  evasion  turns,  are  calculated, 

(See  Figure  3,)  The  time  to  turn  (TTT),  evasion  turn  in  degrees  (AB), 
and  the  turning  radius  (TR)  are  needed  to  compute  X(N)  and  Y(N).  The 
equations  for  TTT,  X(N),  and  Y(N)  are  computed  as  follows; 

_  arc  length  _  TR  x  AB 
sub  speed  U3 

If  TI  >  TB,  then 

X(N)  *  TR  x  SIN(AB)  +  U3  x  (TB-TTT)  COS(AB)  (27) 

Y(N)  *  TR  -  TR  xCOS(AB)-*-  U3  x  (TB-TTT)  SIN(AB)  ,  (28) 

If  TI  <  TB,  then 

X(N)  ~  TRx  SIN(AB)+U3x  (TB-TTT)  COS(AB)  +  (TB-Tl)U2xC0S< AB)  (29) 

Y(N)  -  TR-TRxCOS ( AB ) +  U3  x (TB-TTT) SIN(AB) + (TB-Tl) 

U2  x  SIN(AB) #  (30) 

The  time  Tl  is  the  length  of  time  the  submarine  remains  at  fast  speed  (U3); 
after  which  it  changes  to  silent  speed  (U2).  Equations  (29)  and  (30) 
account  for  the  speed  change  occurring  before  TB  has  terminated.  The 
point  (X(N) , Y(N))  can  be  considered  the  submarine's  actual  position  after 
the  N  evasion  turn  is  taken. 

The  estimated  position  EP  is  the  predicted  position  of  the  submarine 
after  the  interval  of  blind  time  has  elapsed*  Therefore,  EP  is  (DIST,0)* 
The  distance  between  each  (X(N),Y(N))  and  EP  constitutes  a  bias  (BIAS(N))* 
This  is  an  important  factor  in  computing  the  probability  of  success.  The 
bias  for  the  N^1  evasion  turn  is: 

BIAS(K)  =  "Y  (X(N)-DIST)2  ■+-  <Y(N))2  .  (31> 

As  in  the  unalerted  case,  the  parameters  BIAS,  EXPRNG,  and  SIGJT  must  be 
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computed  (EXPRNG  and  SXGJT  are  the  same  as  in  the  unalerted  case)  in 

order  to  use  the  circular  coverage  function.  The  circular  coverage  func 

tion  is  used  to  compute  the  probability  (PPK(N)),  for  each  N,  that  the 

vehicle  can  regain  contact  at  EP  given  the  total  standard  deviation  of 

error  (SIGJT),  the  bias  (BIAS(N))  between  the  terminal  point  of  the 

evasion  and  EP,  and  the  expected  active  sonar  range  (EXPRNG). 

The  two  ratios  of  the  parameters  that  are  needed  for  the  circular 

coverage  function  are  and  .  With  these  ratios  the  proba- 

SIGJT  SIGJT 

bilities,  PPK(N) ,  are  computed  for  each  N,  and  then  all  PPK(N) ,  N  s  1 , 
...,  NU+1,  are  summed  and  averaged  over  all  NU+1.  The  average  proba¬ 
bility  is: 

NU+1 

PPK(N) 

PAVG  =  hTT _  (32) 

NU  +  1 


It  might  be  appropriate  to  use  the  option  to  weight  each  of  the  eva 
sion  turns.  If  It  is  believed  that  the  turns  within  the  set  of  evasion 
turns  are  not  equally  likely,  then  the  probability  vector  LTNU(N),  Nml, 
. . . ,  NU+ 1  can  be  used  which  gives  probabilities  that  each  of  the  N 
evasions  is  actually  taken.  If  they  are  not  equally  likely,  then  the 
"weighted  probability  of  success  (WTPAVG)  against  an  alerted  submarine" 
is  computed.  The  equation  for  WTPAVG  is 


(PPK(N)  x  UNU(N) ) 


(NU-+1) 


(33) 


The  results  of  the  mathematical  development  have  produced  methods 
of  computing  the  probability  of  success  (PAVG)  for  an  alerted  submarine 
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with  equally  likely  evasion  turns,  the  probability  of  success  (WTPAVG) 
for  an  alerted  submarine  with  evasion  turns  that  are  not  equally  likely, 
and  finally  the  probability  of  success  (PROB)  for  an  unalerted  submarine* 
These  probabilities  are  used  to  analyze  and  compare  ASW  vehicles* 
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3*  Input  Data  Rules* 

A  complete  run  for  the  computer  program  (Appendix  C)  consists  of 
25  input  cards.  Each  input  data  card  represents  the  value  of  one  vari¬ 
able  except  for  the  last  two  cards,  one  of  which  contains  the  values  for 
the  AA(1) ,  i“l,  .  5  and  the  other  which  contains  the  UNU(N)  values, 

N-l ,  NU+1  * 

One  of  the  first  functions  of  the  program  is  to  set  the  indices 
NO-1  and  M~l,  (See  Appendix  B--Flow  Charts,)  The  index  NO  is  used  to 
number  the  pages  of  the  program  output.  The  index  M  is  compared  with 
NEXT  (initially  set  equal  to  one).  The  comparison  dictates  whether  the 
program  will  continue  or  terminate.  The  value  for  NEXT  determines  the 
number  of  runs  to  be  executed  and  is  an  input  to  the  program.  The  compu¬ 
ter  then  reads  the  input  data  cards  from  SUBROUTINE  INPUT  after  which 
the  25  values  are  printed  on  page  1  of  the  output.  The  index  NO  is  in¬ 
creased  by  one.  The  programmed  computations  are  performed  and  the 
results  are  printed  on  page  2  of  the  output  (Figure  6),  The  indices  NO 
and  M  are  then  increased  by  one.  The  index  M  is  compared  to  NEXT  and 
if  NEXT  is  greater  than  M,  SUBROUTINE  INPUT  searches  for  more  input  data 
cards.  Since  the  program  is  used  to  compare  and  evaluate  systems,  only 
a  few  variables  will  be  changed  for  each  successive  run.  However,  more 
or  all  the  variables  can  be  changed  as  long  as  the  input  procedure  is 
followed.  For  example,  to  compare  the  three  vehicles  in  Table  1,  first 
place  the  parameters  for  the  High  Speed  Surface  Vehicle,  RAI  s  5.0,  and 
a  value  for  NEXT  in  the  data  input  for  Run  1,  (See  Figure  4*)  The  next 
run  consists  of  a  Locator  Card  which  initiates  the  reading  of  Variable 
Change  Card,  which  in  turn  changes  the  value  of  RAI  to  10,  The  subsequent 
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card  would  be  the  Information  Card>  which  indicates  no  further  changes 
for  the  second  run  --  start  computations.  For  this  run  the  variable  RAI 
has  changed  and  all  other  parameter  values  remain  the  same. 

If  another  variable  change  were  desired ,  another  Locator  Card, 
Variable  Change  Card,  and  Information  Card  would  be  placed  after  the 
RAI  Variable  Change  Card  as  indicated  in  Run  3  in  Figure  A, 

A  run  following  the  initial  run  takes  a  minimum  of  three  cards.  If 
only  one  variable  is  being  changed  (e.g.  RAI  s  10),  then  three  cards  are 
necessary*  This  procedure  is  followed  for  RAI  a  5,  10,  . 60  for  each 
system  in  Table  1  in  order  to  compute  probabilities  of  success  for  each 
vehicle  over  the  entire  RAI  range. 
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TABLE  1 


PARAMETER  INPUTS 
FOR  SAMPLE  ASW  SYSTEMS 


HIGH  SPEED 
SURFACE  VEHICLE 


VERTOL 

VEHICLE 


SEAPLANE 
TYPE  VEHICLE 


ro 


U1 

8.0 

SLO  SUB  16.0 

FAST  SUB 

U2 

3.0 

8.0 

(SAME  SUB  SPEEDS  FOR  THE 

U3 

15.0 

30.0 

Z 

45.00 

120.00 

PTOT 

3.00 

4.0  0 

TOT 

.00 

1.00 

CLOT 

.00 

2.00 

RELT 

3.00 

5.00 

T  1 

10.00 

10. CO 

THETA 

45.00 

45.00 

BETA 

10.00 

10.00 

TR 

40C0.00 

4000.00 

PS  I 

60.00 

60.00 

DELTAL 

330.00 

330.00 

BRGER 

1  .00 

2.00  . 

T 

EXPRNG 

3.00 

5.00 

PRCNTN 

4.00 

2.00 

PRCNTS 

5.00 

5.00 

THDEL 

5.00 

5.00 

UERR 

5.00 

2.  CD 

NU 

6 

6 

~AA  1 1 ) 

-1.28  .  -.52  .  . 

00  ,  .52  ,  1.28 

(SAME 

FOR 

3 

UNUtN) 

.  10 

t  *05  f  *00  t 

.10  ,  .40  ,  .30  ,  .05 

(SAME 

FOR 

3 

3  SYSTEMS) 

20C.C0 
5.00 
5.00 
5. CO 
5.00 
10.00 
45. CO 
1C. CO 
4C0D.00 
60.00 
330. CO 
.50 
6.00 
4.00 
5. CO 
5. CO 
5.00 
6 

SYSTEMS) 

SYSTEMS) 


4*  An  Example  and  Some  Applications, 

The  simulation  possesses  the  capability  o£  being  used  as  a  tool  to 
analyze  and  compare  ASW  vehicles  as  explained  in  the  previous  sections. 

In  order  to  illustrate  this  capability,  the  parameters  listed  for  each 
system  in  Table  1,  and  ranges  RAX* 5 ,  10,  15,  . ..  ,  60  nautical  miles 
were  used  as  input  data  in  the  computer  program  (Appendix  B)  and  pro¬ 
cessed  in  the  CDC-1604  computer.  The  program  is  written  in  FORTRAN- 60. 
The  flow  chart  is  illustrated  in  Appendix  C. 

A  sample  output  for  the  Seaplane  Type  Vehicle  using  the  parameters 
of  Table  1  and  RAI*45  is  shown  in  Figures  5  and  6.  The  former  is  a  list 
of  all  the  parameter  values  for  the  indicated  run.  This  list  enables 
one  to  insure  that  the  correct  parameters  for  that  run  are  being  properly 

-c- 

read  into  the  computer.  The  latter  page  contains  the  probabilities  of 
success  for  one  run.  The  probabilities  for  each  RAX  are  plotted  on 
graphs.  Each  vehicle  is  depicted  on  an  individual  graph.  Each  graph 
has  six  curves  --  three  drawn  from  values  calculated  from  using  a  slow 
submarine  as  an  adversary  and  three  drawn  from  using  a  fast  submarine  as 
an  adversary  (see  Figures  7,  8,  9). 

The  parameter  values  in  this  sample  are  strictly  hypothetical  and 
demonstrate  only  the  use  of  the  simulation.  However,  with  these  chosen 
parameter  values  it  can  be  seen  in  Figures  7,  8,  and  9  that  throughout 
the  range  of  RAI  the  High  Speed  Surface  Vehicle  has  probabilities  of 
success  considerably  lower  than  the  other  systems.  The  slower  the  speed 
of  the  vehicle,  the  lower  the  probability  because  of  a  low  closing  rate 
and  long  blind  time*  The  probabilities  of  success  are  higher  at  shorter 
RAI  for  the  VERTOL  Vehicle  as  compared  with  the  Seaplane-Type  Vehicle, 
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OUTPUT  RUN  5  OF  22 


PAGE  1 


PARAMETER 

VALUE 

UT 

8.00 

U2 

i  .00 

U  3 

15.00 

Z 

200.00 

PTOT 

5  .00 

TOT 

5  .00 

CLOT 

5.00 

RELT 

5. CO 

T  1 

1C.  00 

THETA 

45.00 

BETA 

10.00 

R  A  I 

45.00 

TR 

40C0.00 

PRC  NT  N 

4.00 

PRC  NT  S  . 

5.00 

BRGER 

.50 

EXPRNG 

6.00 

PSI 

60. CO 

DELTAL 

330.00 

NU 

6 

NEXT 

24 

UERR 

5. CO 

THDEL 

5.00 

COMPLETE  INPUT  PARAMETER  SET 


MEANING 

SUBMARINE  CRUISE  SPEED  (KNOTS) 

SUBMARINE  SILENT  SPEED  (KNOTS) 

SUBMARINE  FAST  SPEED  (KNOTS) 

VEHICLE  TRANSIT  SPEED  TO  DATUM  (NEGLECT  ACCEL.  CR  DECEL.) 

PRO  TAKE  OFF  TIME  (MIN) 

TAKE  OFF  TIME  (KIN  I 
LAND  TIME  (MIN) 

TARGET  RELOCATION  TIME  (MIN) 

ELLAPSED  TIME  -OF  SUB  AT  SPEED  U3  (MIN) 

COURSE  OF  SUB  AT  LAST  KNOWN  POSITION  (DEGREES) 

TRUE  BEARING  OF  SUB  FROM  VEHICLE  (DEGREES) 

RANGE  TO  SUP  AT  LAST  KNOWN  POSITION  FROM  VEHICLE  (MILES) 

TURNING  RADIUS  OF  SUB  (VOS) 

NAV.  ERROR  IN  PERCENT  DISTANCE  TC  SUB 
SENSOR  RANGE  ERROR  IN  PERCENT  SENSOR  RANGE 
SENSOR  BEARING  ERROR  IN  DEGREES 

RANGE  EXPECTED  FROM  SENSOR  FOR  LOCALIZATION  AT  DATUM, MILES 

LEFT  BOUND  FOR  SUB  EVASIVE  TURN  IN  DEGREES 
RIGHT  BOUND  FOR  SUB  EVASIVE  TURN  IN  DEGREES 
NUMBER  OF  INCREMENTS  CONSIDERED  IN  PSI-CELTA  RANGE 
NUMBER  OF  SETS  OF  INPUT  DATA 
SUBMARINE  SPEED  ERROR (KNOTS) 

SUBMARINE  COURSE  ERROR  (DEGREES) 


FIGURE  5 


OUTPUT  RUN  5  OF  22 


PAGE  18 


INPUT  PARAMETERS 


TARGET 


VEHICLE 

PRE  TAKE 

TAKE  OFF 

LAND 

RELOCATION 

TIME 

SPEED 

OFF  TIME 

TIME 

TIME 

(KNOTS ) 

(MIN  ) 

(MIN) 

(  M 1  N  ) 

(MIN) 

200 .CO 

5.00 

5.  OC 

5.00 

5.  00 

SURMRI  NE 
CRUISE 
SFEEU 
{ KNCTS ) 

8.00 


SUBMARINE 

SILENT 

SPEEC 

(KNOTS) 

3. CO 


SUBMARINE 
FAST 
SPEER 
(KNOTS ) 

15. OC 


SU  PMARINE 
TO  VEHICLE 
RANGE 
(MILES) 

45.00 


SUBMARINE 

COURSE 

(DEGREES) 


TRUE  REARING 
SUB  FROM 
VEHICLE 
(DEGREES) 


NAVIGATION 
ERROR 
(PERCENT) 
DISTANCE ) 


1*5  .CO 


10.00  '4.00 


SENSCR  EXPECTEC 

RANGE  SENSCR  RANGE  CF 

ERROR  BEARING  LOCALISATION 

(PERCENT)  ERROR  SENSOR 

RANGE)  (DEGREES)  (MILES) 

5.0C  .50  6.00 


N> 

■O 


AVERAGE  SUCCESS  PROBABILITY  FOR  7  EVASION  TACTICS 
AVERAGE  WEIGHTED  SUCCESS  PROF' ABILITY  FOR  7  EVASION  TACTICS 
.  AVERAGE  SUCCESS  PROBABILITY  FOR  UNALERTED  SUBMARINE 


IS 

.539 

IS 

.  C7  5 

IS 

.  679 

COMPUTED  VALUES  FOR  SUBMARINE  EVASION 


- 

BLIND 

DISIANCE  TRANSIT 

RELATIVE 

EVASIVE 

EVASION  TURN 

SUCCESS 

TIME 

ro  CATUM 

SPEED 

SIG  TOTAL 

NUMBER 

( DEGREES ) 

PRCSABIL ITY 

(  M I  N  ) 

(MILES) 

( KNCTS ) 

(MILES) 

1 

-50.00 

.539 

33.96 

52.18 

193.59 

2.03 

c 

-15. OC 

.  621 

33.96 

52.  1  E 

193.39 

2.03 

.00 

.656 

23.96 

52.18 

193.39 

2.03 

4 

15.00 

.  621 

33.96 

52.  18 

193.39 

2 .03 

5 

3C.0C 

.559 

53.96 

52.18 

193.39 

?  2.02 

6 

45.  OC 

.441 

35.96 

l  52.18 

1 93.39 

2.03 

7 

60.00 

.  555 

33.96 

52.18 

192.39  * 

2.03 

t 

FIGURE 

6 
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V  *  no 

RT  s  il 

SR6ER  -i,o 
PRCNTN  =2.0 
EXPRNG  sS.O 
PROMTS  =  S.O 
7HDEL  =5.0 
UERR  =1.0 
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but  a  crossover  point  occurs  at  approximately  RAI=25,  where  the  Seaplane* 
Type  Vehicle  has  higher  probabilities  of  success*  The  crossover  is  a 
result  of  the  values  selected  in  the  sample  for  the  vehicle  speed,  reac¬ 
tion  time,  and  expected  sonar  range*  From  these  results  the  best  vehicle 
can  be  chosen  and  the  other  vehicles  investigated  for  improvements* 

In  addition  to  the  example  stated  above  several  applications  are 
suggested  for  the  simulation*  One  application  is  to  use  the  simulation 
to  determine  what  parameters  or  characteristics  the  vehicle  and  the 
sensors  must  have  to  obtain  a  given  probability  of  success.  For  example, 
if  the  state-of-the-art  for  active  sensors  dictates  that  five  nautical 
miles  is  the  maximum  range  expected  for  active  sensors  for  the  next  ten 
years,  then  it  would  be  advantageous  to  find  values  of  vehicle  speed, 
reaction  time,  navigation  errors,  and  prediction  errors  that  can  be  tol¬ 
erated  for  the  vehicle  and  still  maintain  a  given  probability  of  success* 
The  development  of  a  new  ASW  vehicle  could  be  directed  towards  achieving 
these  parameter  values* 

Another  application  can  be  made  if  a  new  sensor  is  being  proposed 
for  a  given  vehicle*  The  values  for  the  expected  active  sensor  range 
and  the  sensor  errors  can  be  computed  for  a  given  probability  of  success* 
These  parameter  values  can  be  used  as  design  limits  in  the  development 
of  the  proposed  sensor* 

Finally  the  simulation  can  be  used  if  intelligence  indicated  that 
a  new  enemy  submarine  is  in  existence*  The  simulation  can  be  used  to 
indicate  how  well  the  present  vehicles  will  perform  against  the  new 
threat  relative  to  the  old  threat. 
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5. 


Conclusions . 


The  example  in  Section  4  is  the  type  of  comparison  that  can  be  per¬ 
formed  for  various  kinds  of  ASW  vehicles*  There  are  several  decisions-- 
evasion  tactics  used,  submarine  type  used,  and  so  on- -that  have  to  be 
made  prior  to  using  the  program;  but  once  the  decisions  are  made  for  a 
given  situation,  they  have  to  remain  the  same  for  all  systems  being  com¬ 
pared*  Several  modifications  and  suggestions  for  the  model  are  discussed 
below. 

It  might  be  argued  that  an  elliptical  coverage  function  should  be 
used  instead  of  a  circular  coverage  function  because  of  the  elliptical 
sensor  locating  error.  This  change  can  be  made  by  rewriting  SUBROUTINE 
OCIP  of  the  computer  program  (Appendix  C). 

Instead  of  evasion  turns  limited  to  90  degrees  either  side  of  the 
course  line,  they  could  be  programmed  for  360  degrees  evasion;  i*e* ,  180 
degrees  either  side  of  the  course  line* 

It  must  be  realized  that  the  depth  of  the  submarine  has  no  direct 
influence  on  this  simulation*  Therefore,  this  has  to  be  kept  in  mind 
when  parameter  values  such  as  EXPRNG  are  put  into  the  simulation.  It  is 
possible  that  EXPRNG  could  be  less  in  an  evasion  problem  where  the  sub¬ 
marine  is  more  likely  to  increase  its  depth,  and  therefore,  EXPRNG  should 
be  reduced* 

It  would  be  possible  to  give  the  vehicle  freedom  of  movement  before 
transit  commences,  rather  than  being  motionless.  But  the  advantages  of 
the  passive  tactic  for  the  vehicle  are  then  relinquished. 

The  restriction  on  Tl  in  the  alerted  submarine  case  could  be  abol¬ 
ished  by  programming  for  the  possibility  of  having  a  speed  change  exe¬ 
cuted  before  a  turn  is  completed. 
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Possibly  a  passive  sensor  detached  from  the  vehicle  could  transmit 
current  information  about  the  submarine's  actions  while  the  vehicle  is 
inbound  to  intercept  the  submarine*  This  would  increase  the  probability 
of  success* 

In  the  foregoing  sections,  the  simulation  has  been  described  and  the 
mathematical  development  shown.  In  addition,  the  use  of  the  computer 
program  has  been  illustrated  and  an  example  and  its  results  briefly  ana¬ 
lyzed.  The  use  of  the  simulation  is  simple  and  the  results  easily  under¬ 
stood*  Therefore,  it  is  recommended  that  the  use  of  this  simulation  be 
employed  to  compare  and  analyze  proposed  ASW  vehicles  and  operational 
ASW  vehicles. 
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APPENDIX  A 


GLOSSARY 

Notation 

Meaning 

A(l) » * ■ • >A(5) 

*  Number  of  standard  deviations  from  the  mean  of 
the  normal  density  function.  Each  A(i)  is  the 
midpoint  of  a  segment  of  area  which  represents 
20%  of  the  total  area  under  the  normal  density 
function.  The  A(i)’s  are  used  in  the  predic¬ 
tion  error  of  the  unalerted  submarine  case. 

AB 

Evasion  turn  in  degrees  from  the  estimated 
submarine  course. 

BETA 

True  bearing  of  the  submarine  from  the  vehicle 
(degrees) , 

BIAS(N) 

-  Distance  from  EP  to  the  terminal  point  of  the 
evasion. 

BRGER 

Sensor  bearing  error  in  degrees. 

CLDT 

-  Land  time  (min.)  which  includes  time  to  land 
into  the  wind  and  taxi  time  to  get  Into  posi¬ 
tion  to  employ  active  sensors. 

DELTA 

Bound  to  the  right  of  the  known  submarine 
course.  This  bounds  to  the  right  the  set  of 
possible  evasion  turns  (degrees). 

DIST 

Distance  the  submarine  travels  during  blind 
time,  assuming  a  constant  course  and  speed. 

EP 

Estimate  position  of  the  submarine  after  blind 
time  has  expired. 

EXPRNG 

-  Expected  range  of  the  active  sensor  used  by  the 
vehicle  in  re-establishing  contact  at  the  EP 
(in  nautical  miles). 

NEXT 

Number  of  data  runs  for  the  computer  program. 

NU 

-  Number  of  increments  considered  with  the 

(PSI-DELTA) ■ range  of  the  set  of  possible  turns. 

PAVG 

-  Average  probability  of  success  for  N  evasions 
considered  in  the  alerted  submarine  case. 

PRCNTN 

-  Navigation  error  in  percent  of  distance  trav¬ 
eled. 
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Notation 

Meaning 

PRC NTS 

-  Sensor  range  error  in  percent  of  sensor  range* 

PROB 

*  Probability  of  success  for  the  unalerted  sub¬ 
marine  case. 

PSX 

*  Bound  to  the  left  of  the  last  known  submarine 
course*  This  bounds  to  the  left  the  set  of  pos 
sible  evasion  turns  (degrees)* 

PTOT 

Pre  take-off  time  (min.)  which  commences  when 
the  vehicle  can  no  longer  observe  the  submarine 
and  includes  decision  time,  warm-up  time,  and 
sensor  retrieval  time. 

RAI 

-  Range  from  the  vehicle  to  the  submarine  at  the 
last  known  position  (nautical  miles). 

RELT 

-  Target  relocation  time  (rain)  which  includes  the 
time  to  deploy  sensors  and  ends  when  the  sub¬ 
marine  is  relocated. 

RT 

*  Reaction  time  which  is  the  sum  of  PTOT,  TOT, 
CLDT,  and  RELT. 

SBIAS(I,J) 

Distance  between  the  (i,j)^  cell  and  EP* 

SIGJT 

-  Total  standard  deviation  of  errors. 

TB 

-  Blind  time  which  is  the  sum  of  the  reaction 
time  (RT)  and  the  transit  time  (TD) . 

TD 

-  Time  it  takes  the  vehicle  to  travel  from  Its 
initial  position  to  EP. 

Tl 

*  Tactical  interval  of  time  (min*)  after  evasion 
commences  during  which  the  submarine  travels  at 
fast  speed.  After  Tl  has  elapsed,  the  sub¬ 
marine  slows  to  silent  speed. 

THDEL 

-  Submarine  course  error  (degrees). 

THETA 

Course  of  the  submarine  (degrees)  at  the  last 
known  position. 

TOT 

*  Take-off  time  (min.)  which  includes  take-off 
roll  and  the  time  to  maneuver  into  the  wind. 

TR 

-  Turning  radius  of  the  submarine. 
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Notation 

Meaning 

UERR 

Submarine  speed  error  (knots). 

UNU( l) , . - . , UNU(NLH-  l)  - 

The  set  of  probabilities  for  N  =  l ,  *  * . ,  NU+1 
evasion  turns.  These  are  used  when  evasion 
turns  considered  are  not  equally  likely. 

Ul 

Submarine  cruise  speed  (knots). 

U2 

Submarine  silent  speed  (knots). 

U3 

Submarine  fast  speed  (knots). 

V 

Average  vehicle  transit  speed  (knots) 

WTPAVG 

Weighted  probability  of  success  for  the  alerted 
submarine  case. 

(X2,y2) 

Last  known  position  of  the  submarine. 

(X(N).Y(N)) 

Terminal  point  of  the  evasion  after  blind 

time  has  expired. 

2 

Average  vehicle  transit  speed  (knots).  Neglect 
acceleration  and  deceleration  time. 
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APPENDIX  B 


INPUT  AND 
UNIT  CHANGES 


3' 5 


MANEUVERING  BOARD 


36 


37 


38 


ALERTED  CASE 
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40 


i 


41 


J 


I 


42 


UNALERTED  CASE 


43 


SUBROUTINE  INPUT 


-P' 

U1 


PROGRAM  ASWSYS 

C  THE  PROGRAM  CUMPUTLS  PROBABILITIES  OF  AN  ASW  VEHICLE  VS.  AN 

C  EVASIVE  SUBMARINE  AS  A  FUNCTION  OF:  REACTION  OF  THE  VEHICLE. 

Cl  MEN SION  Y3( 2) ,X3(2) , DELTA! 3C  > ,X( 301 ,Y(iO )  , BIAS (301 , AA( > ) , 

ISUCXt 5, 5) t SUOY ( 5i  5 ) , SB  IAS ( 5, 5 ) ,PPK< 1 5 ) , UNU( 15)  ,  ZZ( 4 ) 

NEXT=2 
.  N0=  1 

M=  1 

20  IF  (NEXT-M  )  9CGC-,  30 , 30 

30  CALL  INPUT! Ul« U2,  U3 ,  Z ,  PTOT  , TOT, CLDT, RSLT, T1 , THETA,  RETA,RAI,  TR, 
1PRCNTN,PRCNTS, BRGER, EXPRNG.PSI  ,UERR ,  THL'E  L  ,  NU  ,  NEXT,  CELT  AL,  SUM,  AA,  IF 
2LAG  , UNUi D  EL  N » SUMM ,M ) 

C  OUTPUT  ALL  INPUTS. 

. GO  TO  1999 

C  CONVERT  KNOTS  TO  MILES  PER  MINUTE,  AND  DEGREES  TC  RADIANS 

4C  SV=Z*. 01667 
SU1=U 1*.0 1667 
$U2=U5*. 01667 
SU3=U2* • 0 1 667 
SUERR=UERR*.C1667 
STH£TA=TH ETA/ 57 .29573 
SBCTA=8ETA/S7. 29578 
S  P  S  I  =  P  S I  / 5 7 . 2  9  57S 
SOELT A = CELT AL 757.2957 G 
S8RGER=BRGER/57. 29578 
STHCEL=THUEL/57 .29573 

C  SET  UP  MANEUVERING  BOARD  SOLUTION  TO  SOLVE  FOR  THE  .PCSITCN  THE 

C  SUBMARINE  WILL  BE  WHEN  THE  VEHICLE  INTERCEPTS  IT,  ASSUMING  A 

C  CONSTANT  COURSE/SPEEC  FOR  THE  SUBMARINE. 

Xl=SUl*SINF(STHETA) 

Y1=SU1*CDSF<  STHETA  ) 

C  ( X  I , Y 1 )  AND  (0,0)  ARE  TWO  POINTS  THAT  INDICATE  THE  CUS/SPD  OF  THE 

C  SUBMARINE. 

X2=RA I*SINF(SBETA )  / 

Y2  =  RAI-»C0SF(  SBETA  ) 

C  ( X2 , Y2 )  IS  THE  LAST  KNOWN  POSITION  OF  THE  SUBMARINE  RELATIVE  TO 

C  VEHICLE  AT  (C,C).  NOW  CHECK  SPECIAL  CASES. 

C  THE  FOLLOWING  13  STATEMENTS  ARE  USED  TO  ROUND  OFF  Xl,X2,Yl,Y2  TO 

C  ZERO  IF  THEY  ARE  WITHIN  4  CR  -  .001  'OF  ZERO. 

ZZ(  1 ) =X1 
ZZ  (  2 ) =Y  1 
ZZ ( 3 1 =X2 
ZZ  (  4 ) =Y 2 
DC  52  1=1 ,4 

AP.ZZ=ABSF(ZZ(I)  )  .  ■ 

IF(  ABZZ-. 001) 51,51,52 

51  ZZ (  I) =C  . 

52  CONTINUE 


CC10 
CO  20 
0050 
C04C 
0050 
0060 
0070 
CCBC 
C090 
0100 
C  1  10 
C  120 
0130 
0140 
C  1  50 
0160 
0170 
0180 
C  1  90 
0200 
C2  10 
C22C 
C250 
C240 
0250 
02  6  C 
02  70 
02S0 
C  290 
0300 
03  10 
0320 
C3  3C 
0340 
0350 
036C 
C  37  0 
0380 
0390 
C4CC 
C4  10 
0420 
0430 
CULC 
04  50 
0460 
0470 
0430 
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APPENDIX  C 


\ 

\ 


.X  .  _J 


O' 


102 

103 

104 

105 

106 


xi=zzm 

Y1 =ZZ (2  ) 

X2  =  ZZ  ( 5 ) 

Y2=ZZ (4) 

IF(X2 1105, 100, 105 
100  Y3( 1  1  =SCRTF( SV«*2-X1*«2) 

Y3 ( 2 1 =— SQRTF ( SV**  2— Xl **2 1 
Y3A  =  Y  2(1) 

Y3E=Y  3(2) 

AB 1 YY  =  ABS  F  ( Y2-Y3A 1 
AB2YY=ABSF(Y2-Y3B) 

IF (AB  1YY-AB2YY) 102,9050, 103 

GO  TO  104 
1  =  2 

X3 (  I ) =X 1 
GO  TO  2CC 
IF(Y2)110»1 06 ,110 
X3 C  1  )=SQRTF(SV**2-Y1**2) 

X3( 2) =-X3(  1 ) 

X3 A=X  2(1 ) 

X3B=X  3(2) 

AB  1  XX  =  ABSF(  X2-X3A  ) 

AB2XX=ABSF ( X2-X3B ) 

IF( AB  IXX-AB2XX) 107,9050, 1 08 
1  =  1 

GO  TO  109 

I=?  2 

Y3 ( I) =Y 1 

GO  TC  200  J 

B=( (-2.*X1*Y2*Y2)/(X2**2) )+ ( (2 .* Yl * Y2 S /X2 1 
A= 1 .+ (Y2**2/X2«*2 ) 

C=( (Y2«*2*Xl**2)/X2**2)-(2.*Yl*Y2»Xl 1 /X2+Y1**2-SV*«2 
X3 ( 1 ) =( -B+SQRIF (3  **2-4«*A*C ) ) /(2.*A) 

X3( 2  j  =( -B -SQRTF (B**2-4.*A*C) )/(2.*A) 

X3A=X2 ( 1 1 
X3B=X2( 2) 

ABSXXl=ABSr (X2-X3 A ) 

ABSXX2=ABSF(X2-X3B) 

IF ( ABSXX1-ABSXX2)  125,9050,155 
125  1=1 

Y31=SCRTF(SV**2-X3A**2) 

Y32=-Y31 

ABSYY  1  =  ABSF( Y2-Y31 1 
ABSYY2=ABSF(Y2-Y32 ) 

IF ( ABSYY 1-ABSYY2) 130,9050,135 
130  1=1 

Y3(I)=Y31 


107 

108 

109 

110 


0490 
0500 
0510 
C520 
C530 
0540 
0550 
0560 
0570 
0580 
C59C 
0600 
0610 
0620 
C63C 
0640 
0650 
0660 
C670 
0680 
0690 
C7C0 
0710 
0720 
0750 
C740 
07  50 
0760 
C77C 
0780 
0790 
0800 
0810 
0820 
0830 
0840 
C850 
0860 
0870 
0880 
C890 
0900 
0910 
092C 
09  30 
0940 
09  50 
096C 


noon  n  nn 


-j 


.  GO  TO  200 
135  1=1 

Y3(  I ) =Y  32 
GO  TO  200 
155  1=2 

Y33=SQRTFI$V**2-X3G**2) 

Y34=- Y35 

ABSYY1  =  AB$F(Y2-Y33  ) 

ABSYY  2=  AESF (  Y2-Y34  ) 

IF  < A8SYY1-A9SYY2) 14C, 9050,  150 
11)0  1  =  2 

Y3 (I) =Y  33 
GO  TO  200 
150  1=2 

Y3( I ) =  Y34 

200  X3Il=X3t II- 
Y31  I=Y3(  I  1 

RELSPC=SCRTF(  ( X 31  I -X 1 )« *2  +  ( Y 5 1  I - Y 1 1 ** 2 ) 

RSPC=FELSPD/.Cl667 
C  COMPUTE  FLYING  TIME  TO  DATUM 

TD  =  P>AI/RCLSPrj 

C  COM PU T C  REACTION  TIME 

RT  =PTGT+TCT+CLCT+RELT 
IX  =  0 

'  SU  BT  =  R  A I  /  SU  3 
IF l SUeT-RT)202, 202,201 
202  IX=1 

201  CONTINUE 

C  COMPUTE  OVER  ALL  REACTION  TIME  (CALLED  BLIND  TIME). 

TB=TD+RT 

C  COMPUTE  DISTANCE  TRAVELED  BY  SUBMARINE  DURING  RLIND  TIME. 

DIST=SU3*TB 

(XL ,Y4)=PREDICTtO  DATUM  POSITION  OF  SUBMARINE  RELATIVE  TO  VEHICLE 
AT  (0,0.  / 

X4=  D  I  ST  *S  IN  F  (  ST  HE  T  A  )  +  X2 
Y4  =  DI ST  »COS  F ( STHET  A )  +  Y2 
DISTANCE  FLOWN  TO  DATUM. 

DI  SUAT  =  SQRTF{  XL**  2. +Y4  **2  ) 

COMPUTE  PROBABILITIES  FOR  ALERTED  SUBMARINE  EVADING.  PLACE  SLE  IN 
(0,0)  POSITION  HEADING  000.  RANGE  OF  SUB  EVASION  IS  FROM  DELTA 
CEGREES  ro  RIGHT  OF  GOO  TO  »SI  DEGREES  TO  THE  LEFT  WITH  INCREMENTS 
OF  ETA  DEGREES.  FOR  COMPUTING  PURPOSES  CGC  IS  LAST  KNCWN  SUB  CUS. 
RNU=NU 

I F (  1 . S70S-SDELTAJ  204 ,9050, 21 C 
IF  [  1. fTCS-SPSI  )21 0,22  5,205 
FT  A=( SPS I  +  ( 6. 2B31  9-SCtLTA} ) /RNL 
GO  TO  2  12 

ET  A  =  ( SPSI -SDE  LT  A ) / RNU 


2  04 

?C5 

2  1C 


C97C 
C  9  f)  0 
0990 
100C 

10  1C 
1020 
1030 
1  0  4  C 
1050 
1060 
1070 
1080 
1090 
1  100 
1110 
1  120 
1130 
1140 
1  150 
1  160 

1 1  70 
1  1  80 
1190 
1200 

12  10 
122C 
12  30 
12  DO 
1250 
126C 
1270 
1230 
1290 
1300 
1310 
1320 
1330 
1340 
1350 
1  3o0 
1370 
1  380 
1390 
14  CC 
14  10 
1420 
1430 
1440 


212  STR=T K/20CC . 

SAVE  1 =  STR 
COMPUTE  ERRORS 
S  I  GNAV=  .  0  1  *PRCNT\'*RA  1 
S I  GSR=  .  CURRENTS*  RA1 
S I GSh  R=RA I  *S I KF (SBRGER ) 

SI  GAPR=S()RTFt  5 1  GSR  *  S 1 GSI3R  > 

$IGJV=  SCRTF(  SIG'J  AV**2+SI  GAPR**2  ) 
S AVE=TB 

IF (TB-T1 )299,  299,  405 
405  I F ( T 1  >9050,  299,  407 
407  I FL  AG= 1 

TB=T  1 

299  UN=NU+1 
NUN=U  A 

00  50C  N»1,NUN 
IF (N- 1>9050, 301 ,3 10 
301  C£LN=SI)£LTA 
210  DN=AOSF (CELN-d .2832 1 

IF ( DN-. 0009 >312,3 12, 314 

312  C£LN=  6 . 2832 

314  lF(0ELN-6.283w)300,31 5, 666 

315  CONTINUE 
CELTAIN )  =  0. 

TTT  =0  . 

AB  =  0  . 

ST  R=0 . 

BB  “0  . 

GO  TO  35C 

666  DELN=CELi\l-6.2B32 

300  DM=ABSF(OELN- i .5703) 

IF ( CM-. 0009)3 13,312, 316 

313  DELN=  1 . 57C3 

316  IF  (  D'ELN-1  .  570  3)320,230,540 
320  TTT=STR*0ELN/SU3 

AD  =DE  LN 

DELTA  IN >=DELN*57.2957S 
B  B  =  A  B 
GO  TO  350 
330  AB= 1.5708 

DELTA ( N ) =90 . 

X( N)=S1R 

YIN  )  =  STR+  t  TE-TTT) *SU3 
GO  TO  1^09 

340  AB=- ( 6 . 28  32-DEL N) 

DELTA (N)=-t  6. 2852-DFLN) *37.29  578 

TTT=STR*AB/SU3 

RB=  AB 


1450 
1460 
1470 
1480 
1490 
1500 
15  1C 
1520 
1530 
1  56  C 

15  50 
1560 
157C 
1530 
1590 
160C 

16  1C 
1620 
1630 
1640 
1650 
1660 
1  6  70 
1  68  C 
169C 
1700 
1710 
1  72  C 

17  20 
1740 
1  7  5  C 
1760 
1770 
1780 
179C 
1800 
1810 
1820 

18  30 
1840 

18  50 

ieot 

1870 
1880 
1  89C 
1  9CC 

19  10 
1920 


non  no  o  no 


STR=-STR 
TTT=-TTT 

350  X(N  )=STR*  SINF  (  nB)+SU3*  <  TD-TTT  )*C0SF ( Ah ) 

Y(  N)=STR-STR*CCSF  (  BS ) + SU3* ( TB-T T T ) * S I NF < AB  ) 

4  09  I  F  (  IFLAG)  9C5'J,U00, 4  10 
4  10  X  X  X  =  X  (  N  I 
YYY=Y (N) 

T3=SAVE-T 1 

X( N)=TO«SU2*COSF( AR) 

Y(  N  )=TB*SU2*SINF(  AB  ) 

X( N)=XXX+X(N) 

Y(N)  =  YYY+Y(N  > 

40C  DELN=OELN+ETA  • 

XX=0I ST 
YY  “0 . 

COMPUTE  BIAS  FOR  EACH  (X(N),Y(N)}  RELATIVE 
XNX=X (N) 

YNY  =Y  (N  ) 

B I  AS  I N  )  =  SUR  T  F ( (XNX-XX >*«2+ (YNY-YY) **2  I 

D=BI AS ( N ) /SICJT 

AKILL  =  EXPRNG/SIGJT 

STR=SAVS1 

TC=T  j 

COMPUTE  PROBABILITY  FOR  EACH  EVASION  BY  SOLVING 
FUNCT ION. 

CALL  GC IP [ 0, AKILL , PX > 

PPKt  N  )=PK 

SO  MM  =  SUMM  +  PPK  (  N)«-UNU  IN  ) 

500  SUM=SUM+PPK( N ) 

COMPUTE  PROBABILITY  FOR  EVASIVE  SUBMARINE  USING 
UT  P  AVG=SUMM/Ufi 
P AVG- SUM/ UN 
COMPUTE  PROBABILITY 
PRESENCE. 

SU  MUNU=0 . C 
CO  50 2  N= 1 , NUN 
502  SUM UNL=  SU KUNU ■* UN!J  (  N  ) 

LAST  KNOWN  POSITION 
ME  POSITION  AFTER  BLIND 
POSITIONS  OF  CENTERS  Cf 
SU  M=0  . 

00  60C  1=1 
T  B  =  S  A  V  E 
VV  =  (SU-AA(  I  ) 

00  60 C  J=  T  t  5 
THLT  E  R  =  AA I J ) * 

SUBXt 
SU  BY  ( 


TO  (XX, YY) 


CIRCULAR  CCVERAGE 


WE  I GHT ED 


19  30 
19LU 
T95C 
1960 
1970 

19  50 

me 

2CCC 

20  1  0 
2020 
2C3C 
20  40 
2050 

206  C 

207  0 
2030 

20  90 
2  1  C  C 

21  10 
2120 
2130 
2  140 
2150 
2  1.60 
2170 
2130 
2190 
2200 
22  1C 

PROBAR  IL  IT  I  ES  22  20 
22  50 


ASSUMING  SUBMARINE  NUT  ALERTED  GF  VEHICLES 


CF 


SUBMARINE  IS  THE  ORIGIN. 
TINE  (WITH  NO  EVASION) 
CELLS. 


FR EDICT  ED 
IS  (XX , YY I 


SU  E  M  AR  I 
COMPUTE 


,5 


jSUERR  )  «T5 


THDEL 

I,  j  )=VV*C0SF(  THE  TER  ) 
I  , J ) =VV*  SINF ( THETER) 


2?  4  C 
2250 
2260 
22  70 

22  2  C 
2290 
2300 

23  10 
2320 
2  3  30 
2  3^0 
233  C 
2360 

23  70 
2  3  30 
2590 

24  0  0 


non 


XSUB=SUBX (  I ,  J  ) 

YSUB=  SUBY ( 1 ,J ) 

SB  I  AS (I,J>  =  SQRTF( { XX-XSUB)«*2+(YY-YSUB)**£) 

C-SBl  AS!  I  ,  JJ/SIGJT 
AKILL=EXPRNG/SIGJT 
CALL  CCIP  ( Of  AK ILL »  PK ) 

60C  SUM=SUM+PK 
PR08=SUM«.04 
GO  1C  501 C 

PRO  B  IS  PROS  At)  ILI T  Y  THAT  VEHICLE  RELOCATES  SLR  AT  LATUM  IF  SLR 
COES  NOT  EVADE 

PRINT  OUT  INPUT  PARAMETER  SET 
.1999  WRITE  OUTPUT  I  APE  3 , ? 00 1 , M » ME XT , NO 
2001  FORMAT!  lHl, //,5X,  1CHCUTPUT  RUN  , I  3, 1 X ,2HGF , 13 ,79 X, 5HPAGE  ,13) 

WRITE  OUTPUT  TAPE  5,2005 

2005  FORMAT! 1HG.U6X, 28HC0KPLETE  INPUT  PARAMETER  SET) 

WRITE  CUT  PUT  TAPE  3,2010 

23  1C  FORMAT ( ///,20X, 9H  P  A  RAM  E  T  E  R  *  1 IX, 5H VALUE , 15X , 7HMEANI NG) 

WRITE  OUTPUT  TAPE  2,2015, U1 

2015  FORMAT !///.,  2GX ,  2HU 1  ,  1  8X»  F  1C. 4 , 1  OX, 30HSU6MAR INE  CRUISE  SPEEC  (KNOTS 

ID 

WRITE  OUTPUT  TAPE  3,2020, U2 

.2020  FORMAT!/  , 2CX , 2HU2 , 1 8X, F 1 C  .4 ,  1CX , 30HS LBMAR I NE  SILENT  SPEED  (KNOTS) 

u,  1)  _  .  ■ 

o  WRITE  OUTPUT  TAPE  ;t2025<,Ui 

2025  FORMAT!/  ,20X,2HU3,  1 3X , F  1  0 .4  ,  10X, 28 HSUB MARINE  FAST  SPEED  (KNOTS)  ) 
WRITE  OUTPUT  I  APE  3,2030,2 

2030  FORMAT!/  , 20X , 1  Hi ,  1 9X, F I C . 4, 1  OX, 57HVEHIC LE  TRANSIT  SPEED  TC  DATUM 
KNEGLECT  ACCEL.  OR  DECEL  .  )  ) 

WRITE  OUTPUT  1  A  PE  5,2035, PTOT 

2055  FORMAT  (  /  ,  20X  ,  A  HP  TOT  ,  1  6X,  FI  0. 11 ,  10X  ,  23HPRL  TAKE  OFF  TIME  (MIN)) 
WRITE  OUTPUT  TAPE  5,2G4C,TDT 

2040  FORMAT!/  , 20X , 3HT0 T , 1 7X , F 1 0 . 4 , 1  OX , 1 9HT AK£  CFF  TIME  (MIN)) 

WRITE  OUTPUT  1APE  5,^045, CLDT 

2045  FOR  MAT ( /  ,  20X  ,  4HCLCT  ,  16X,F10.4,10X, 15 H LAND  TIME  (MIN)) 

WRITE  OUTPUT  1APE  3,2050, RELT 

^  2050  FORMAT!/  , 20X , 4HR EL T , 16X , F 1 0 . 4 t 1  OX , 2 8HT ARGET  RELCCA1 ION  TIME  (MIN) 

1  ) 

WRITE  OUTPUT  TAPE  3,2C55,T1 

2055  FORMAT ( /  , 2QX , 2HT 1 , 1 8X , F 1 C . 4 , 1 CX , 58HE LEAP SED  TIME  CF  SUB  AT  SPEEC 
1U3  (MIN)) 

WRITE  OUTPUT  I  APE  3, 2060, THETA 

2060  FORMAT!/  , 20X, 5HTHETA, I 5X , FI  0.4 , 1 CX , 46HC0URS E  CF  SUB  AT  LAST  KNCWN 
1  POSITION  (DEGREES)) 

WRITE  OUTPUT  I APE  3,2065, BETA 

2065  FORMAT!/  , 20X , 4HRE TA , 1 6X , FI  0 .4 , 1  OX , 42HTRUE  BEARING  CF  SUB  F  R  CM  VEH 
1ICLE  {DEGREES ) )  ■ 

WRITE  OUTPUT  TAPE  2, 2070, RA I 


24  1C 
2420 
2  4  30 
24  4  0 
245C 
24  60 
247C 
248C 
2490 
2  50C 
2510 
2520 

‘  2530 
2540 
255  C 
2560 
2570 
2580 
259C 
2600 
2610 
262C 
2630 
2640 
2  650 
266C 
2670 
2680 
2690 
2700 
2710 
2720 
2730 
2740 

27  50 
2760 
2770 
2730 

■2790 

2800 

28  10 
28  20 
2830 
2840 
28  30 

2e60 

2870 
2  88C 


J 


207 C  FORMAT  (  /  ,  20X  ,  3  HR  A  I  ,  1  7X  ,  F  1  C  .  4 ,  1  OX,  56HRANGE  TO  SUB  AT  L  AST  KNOWN  PO 
1SITI0N  FROM  VlHICLE  (MILES) ) 

WRITE  OUTPUT  TAPE  3,2080, TR 

2080  FORMAT! /  , 20X , 2HTR , 1 8X, Fl 0. 4 , 1 CX , 27HTURNI NG  RADIUS  CF  SUB  (YDS)) 
WRITE  OUTPUT  TAPE  3 , 2085 » PRC NT N 

2085  FORMAT!/  ,  20X ,  6HPRCNTN ,  1  AX ,  F  I C.  4  ,  1  OX  ,  37HN  A  V  .  ERRCR  IN  PERCENT  CIST 
1ANCE  TO  SUE) 

WRITE  OUTPUT  TAPE  3, 2C90, PRCN TS 

209C  FORMAT!/  » 20X, 4HPRCNTS » I  LX, F 1 0. 4 , 1  OX , 42HS E NSCR  RANGE  ERROR  I  A  PERC 
TENT  SENSOR  RANGE) 

WRITE  OUTPUT  I  APE  3 , 2C95 , BRGER 

2095  FORMAT ! /  , 20X , 5HBRGER , 1 5X , F 1 C  .4 , 10X , 3 1 H SENSOR  BEARING  ERROR  IN  CEG 
1REES) 

WRITE  OUTPUT  TAPE  3 , 2 1 GO , EXPRNG 

2100  FORMAT!/  ,  20X ,  6HE  XPRNG,  1 4X,  F 1 0 . 4 , 1  OX  ,  58HRANGE  EXPECTED  FROM  SENSE! 
TR  FOR  LOCALIZATION  AT  DATUM, FILES) 

WRITE  OUTPUT  TAPE  2,2105, PS  I 

2105  FORMAT!/  ,  20X,  3H?  S  l ,  I  ?X,  F  1C  .  A  ,  10X  ,  42HLEF1  BOUND  .FOR  SUB  EVASIVE  TU 
TRN  IN  DEGREES) 

WRITE  OUTPUT  TAPE  2,2110,DELTAL 

2110  FORMAT!/  , 2CX .6HDELTAL , 1 4X, F 1 0. 4 , 1  OX , 43HR I GFT  BOUND  FOR  SUB  EVAS  IV 
IE  T UR A  IN  DEGREES) 

WRITE  OUTPUT  TAPE  3,2115,NU 

2115  FORMAT!/  , 20X , 2HMU , 1 BX , I  5 , 1 5 X , 47HNUMBE R  INCREMENTS  CCNSIOEREC  IN  P 
1SI-D6LTA  RANGE ) 

WRITE  OUTPUT  TAPE  3,2120, NEXT 

2120  FORMAT  t /  ,20X,4HNEXT , UX, IS, 1 5X , 28HNUMBER  CF  SETS  CF  INPUT  DATA) 
WRITE  OUTPUT  TAPE  3,2125, UERR 

2125  FORMAT!/  ,  PCX  , 4 HU ERR , 16X,F10. 4, 1GX»29HSUH MARINE  SPEED  ERROR  (KNOTS 
T  )  ) 

WRITE  OUTPUT  TAPE  3,2120,THDlL 

213C  FORMAT!/  , 20X , 5HTH0E L , 1 5 X , F 1 C . 4 , T OX , 52HSUBMAR INE  COURSE  ERROR  (DEG 
1  REES )  ) 

1  )  ) 

NO=NO-*  1 
GO  TO  4  0 

9050  WRITE  OUTPUT  TAPE  3, 9055 , K ,NE XT, NO 

9055  FORMAT ( lhl ,/// ,5X,  1 1H0UTPUT  RUN  , I  5 , 1  X , 2H0F , I  3 ,79X , 5FP AGE  ,13,//// 
1 ,5X,i7H£RR0R.  INCORRECT  SOLUTION  BEING  CALCULATED.  CONTINUE  WITH  N 
2 EXT  RLN.) 

NO=NtH  1 
M=M+  1 
GC  TO  2C 

30 TC  WRITE  OUTPUT  TAPE  5, 3C 1 5 , M , N E X T , MO 

3'J  1  5  FORMAT!  1H1  ,//,5X,  1CH0UTPLT  RUN  ,  I  3 , 1  X  ,2H0F  ,  Ii  ,79X,  5HPAGE  ,13) 

WRITE  OUTPUT  TAPE  3,3C20 

3C20  FORMAT! /// ,5GX, 16H INPUT  PARAMETERS ) 

WRITE  OUTPUT  TAPE  3 , 5025 , Z , P TOT , TOT , CEDT , RE L T , U 1 ,U 2 , 1 3 , R A  I ,T NET  A 


28  90 

29  OC 

2  94C 
2950 
2960 
2970 
2930 
2990 

30  GO 
30  tC 
302C 
3030 

3  CHC 
3C50 

30  60 
3070 
3030 
3090 
3100 

31  1C 
3120 
3130 
314C 
3T5C 
3160 
31  70 

31  EC 
3190 
3200 
3210 
322C 

32  30 
3240 
3250 
3260 
3270 
3230 
3290 
3300 
3310 
35.2  0 
336  C 
3  340 
3350 
3560 
3370 
3380 
3390 


U1 

w 


FORMAT  {///,46X,56HTARGET 

1  UP  MAR  I, \E/  IN  ,  :*X,  1  1CH VEHICLE 

2  ION  CRUISE  SILENT 

3H  , 6X , 106HSP£cD  OFF  TIME 
4  C  SPEED  SPEED 

5TS )  (PIN)  (MIN) 

6S )  (KNOTS)  (MILES) 


3025  FORMAT (///, 46X, 56HTARGET  SUBMARINE  SUBMARINE  SUBMARINE  S 

PRc  TAKE  TAKE  CFF  LAND  RELCCAI 

FAST  TO  VEHICLE  SUDMARINE/1 

TIME  TIME  TIME  SFEE 

RANGE  CCURSE/1F  ,5X,110H(KN0 

(MIN)  (MIN)  ( KNOTS)  (KNCT 

(DEGREES) //1H  ,  6X ,F 6 . 2 , 2X , F6 . 2 , 6X , F 
76.2,3X,F6»2»6X,  F6  «2,oX,F6««i,7X,F6.^,6X,F6.^;,oX,F6.2,7X,F6»2) 

WRITE  OUTPUT  TAPE  2 , 5070, BETA , PRCN TN , PRCNTS , BFGER , EXPRNG 
3030  FORMAT ( // / *  49X,  52HSENSCR  EX PECTED / 1 H  ,20X,61HTRUE 

1  BEARING  NAVIGATION  RANGE  SENSOR  RANGE  OF/  1 H  ,21X, 

263HSUP  FROM.  ERROR  ERROR  GEARING  LQCALI2A1I CN 

3/1H  ,23X,59HVLHICLE  (PERCENT  (PERCENT  ERR CR  S 

4ENS0R/1H  , 21X, 59H( DEGREES  )  DISTANCE)  RANGE)  (DEGREES) 

5  ( M I L  E  S ) / / 1 H  »22X,F6.2f9X,F6.2,6X,F5.2,7X,F6.2,6X,F6.2) 

IF ( 1X19050, 3034,3323 
3333  WRITE  OUTPUT  JAPE  3,^C32 

3032  FORMAT!///, 10X, 105H«*  NOTE-RANGE  TO  SUB  IS  CLOSE  ENOUGH  FDR  SUE  TO 
1  CLOSE  VEHICLE  AND  FIRE  WEA-P.CNS  BEFORE  VEHICLE  CAN  ESCAPE.) 

3034  CONTINUE 

WRITE  OUTPUT  TAPE  3, 3 035, NUN , PA VG 

3C35  FORMAT ( ///, 30  <, 32HAVCRAGE  SUCCESS  PROBABILITY  FOR  ,I2,20H  EVASION 
1  T  ACT  ICS  IS  ,  T  .  3 ) 

IF ( SUMUNU ) 9C5C , 30  36, 3C5 1 

3036  WRITE  OUTPUT  TAPE  3,3037 

3037  FORMAT ( / / , 50X , 4SH A VERAGE  WEIGHTED  SUCCESS  PROBABILITY  NOT  USED) 

GO  TO  3039 

30  31  WRITE  OUTPUT  TAPE  3,  30  39 ,  NUN  ,  W  TP  AVG 

3038  FORMAT ( // , 21X, 4 IH AVERAGE  WEIGHTED  SUCCESS  PROBABILITY  FOR  , I2,20H 
1  EVASION  TACTICS  IS  ,F5. 3) 

3039  WRITE  OUTPUT  TAPE  3,3041, PROC 

3041  FORMAT ( /// / , 4  2X  » 5  9HCCMPUT  ED  VALUES  F  CR  UNALERTED  SUBMARINE// 1 H  ,30 
1 X, 28HAVER AGE  SUCCESS  PROBABILITY  , F5.3) 

WRITE  OUTPUT  TAPE  3,3040 

3040  FORMAT ( ///,42X,  57HC0MPUTED  VALUES  FOR  SUBMARINE  EVASION/ If  1 H  ,62 

IX, 3  7H  EL  IN  D  DISTANCE  TRANSIT  RELATIVE /lH  ,17/  ,  9  4  H  E  V  A  S  ION 

2 EVA SI CN  TURN  SUCCESS  TIME  TC  DATUM  SPL 

3CD  SIG  T01AL/1H  , 1 7 X , 94HNUMBER  (DEGREES)  PROEABILI 

4TY  (MIN)  (MILES)  (KNOTS)  (MILES)/) 

DO  30  55-  N  =  1  , NUN 

WRITE  OUTPUT  TAPE  2 , 3050 , N, DELTA  IN) , PPK ( N ) , TE  ,Cl SDAT , RSPC , S I G JT 
3050  FORMAT  (  /  1  H  ,  1  CX ,  I  5 , 3X ,  F8 . 2 ,  1 CX  ,  F 5 . 3 ,  BX ,  Fo  .  2  , £ X  ,  F 6 . 2  ,  11  X  ,  F6 . 2 , 6X ,  F:> 

1.2) 

3055  CONTINUE 
N0=f40+  1 
M=M+ 1 
GO  TO  20 
9000  END 


34  CC 
24  10 
3420 
3430 
34  4  C 
3450 
2460 
3470 
34  90 
3485 
349C 
3500 
2510 
251  5 
3520 
3530 
2540 
3550 
3560 
3570 
3  SCO 
3590 
36CC 
2610 
3620 
3630 
364  0 
2650 
3660 

36  70 
268  C 
3690 
3700 

37  10 
27  20 
2730 
3740 
2750 
3760 
3  7  70 
3780 
3790 
3800 

38  10 
3820 
3830 

384  G 

385  C 


O' 

OJ 


SUBROUTINE  INPUT!  U 1 , U2 , U3 , Z , PTOT , TOT , CLD T , RE L T , T 1 , THE  T  A , E£7  A , R  A  I  , 

1  TR,PRCNTN,PRCNTS,8RG£R ,EXPRNG,PSI ,UERR,THUEL , NU , Nb XT , DELT AL , SLM 
2, A A,  I  FLAG, UNO ,DELN, SLMU,M> 

DIMENSION  AA  (  6 ) ,UNU  (  «  5 ) 


SUM=0 

* 

SUM'M* 

w  * 

IF  LAG 

=c 

DELN- 

c. 

IF  (M- 

1)111,91,93 

91  READ 

1 10, U1 

READ 

1 1C,U2 

READ 

110, U3 

READ 

110,2 

READ 

1 10, PTOT 

READ 

110, TOT 

READ 

110, CLDT 

READ 

11 0, RELT 

READ 

110, T1 

READ 

110, THETA 

READ 

110, BETA 

READ 

1  1  0 ,  RA  I 

READ 

11 0 »  TR 

READ 

1 10, PS  I 

READ 

110, DELTAL 

RL  AD 

1 10 , BRGER 

READ 

110, LXPRNG 

READ 

110,  PRCNl'K 

RE-AD 

no,  PROMTS 

READ 

110, THDEL 

READ 

1  10.UERR 

110  FORMAT!  FI  2 -it  ) 


READ  112, NU 
READ  112, NEXT 
'  1 12  FORMAT ( I  4 ) 

READ  IIS, ( A At  11,1=1,5) 

113  FORMAT! bF 12. 4 ) 

NUN=NU+  1 

READ  1  14,  (UNU(M,N=1  ,NUN  J 
1U!  FORMAT!  12F6.li  ) 

GO  TO  111 
95  REAC.  9l;,  I JK 
-  94  FORMAT! 14 ) 

GO  TOl  1,2,3,  4,S,  6,  7,  8, 9,  10,11 ,12,  15,  14, 15,  16,  17, 13, 19,20,21 ,22, 23‘, 
1  24, 2b ,26,27) ,  IJK 

1  READ  110.U1 
GO  TO  95 

2  READ  1 T 0, U2 
GO  TO  93 


3860 

3870 


3  9CC 

39  1C 
3920 
3930 
3  94  C 
3950 
3960 

3  97  C 
3980 
3990 

40  OC 
4C  10 
4020 
4030 
4040 
4050 
40  AO 
40  80 
4C9C 

4  100 
4110 
4120- 
4  13C 
4  140 
4  150 
4  16C 
4  170 
4  180 
4190 
42  OC 

42  10 
4220 
4230 
4240 
4250 
4260 
4270 
4280 
4290 
4300 

43  1C 
4320 
4  3  30 
4340 


3  READ  11  0, U3 
Gfl  TO  93 

Is  RE  AC  IU),Z 
GO  TO  93 

5  READ  no,PTOT 
GO  TO  93 

6  RE AC  110, TOT 
GO  TO  93 

7  READ  11C, CLOT 
GO  TO  93 

8  REAR  11 0, REL1 
GO  TO  93 

9  RE  AO  11 C , T  1 
GO  TO  93 

1C  READ  110, THETA 
GO  TO  93 

11  READ  1  1  0  f  BET  A 
GO  TO  93 

12  READ  lTO.RAI 
GO  TO  93 

THIS  SPACE  FREE  FOR  AN  INPUT  VARIABLE 

13  GO  TO  93 

1H  READ  1 1 C , T  R 
GO  TO  93 

IE  RE  AC  110, PS  I 
GO  TO  93 

16  READ  I  10, DELTAL 
GO  TO  93 

17  REAR  110, 6RGER 
GO  TO  93 

IE'  READ  110,  EXPRNG 
GO  TO  93 

19  READ  110, PRCNTN 
GO  TO  93 

20  READ  110, PROMTS 
GO  TO  93 

21  READ  110, THDfcL 
GC  TO  93 

22  RE  AC  HO.UERR 
GO  TO  93 

23  READ  11 2, NU 
MJN-NU+  1 

GO  TO  92 

24  READ  112, NEXT 
GO  TO  95 

25  READ  113, (AAI 11,1-1,5) 

GO  TO  93 

.26  READ  11U,  (UNU(M,N  =  1  ,NU.N) 


4350 
4  360 
4  3  7  c 
4  380 
A  390 
4400 
44  1  C 
4  4  20 
4430 
4  44  C 
4450 
4460 
44  70 
44  OC 
4490 
4  500 
4510 
4520 
452C 
4540 
4  5  5  0 
4560 
4570 
4580 
4590 
46  CO 
4610 
4620 

46  2  0 
4640 
4650 
466C 
4670 
4  6  80 
4690 
4  7-CC 

47  3  0 
4720 
4  7  5  C 
4740 
47  50 
476C 
4  7  7  C 
47  3  0 
4790 
4800 
4810 
4820 


i 


s 


w< 

Oi 


GO  TO  93 
27  GO  TO  1 1 1 
111  RETURN 
END 

SUBROUTINE  OC  IP ( A , V , PAV ) 

IFIV-A-5.Q1 1,  1,2 

2  PAV=1.C 
RETURN 

1  IF  ( V-A+5.0  )3,  4,  *4 

3  PAV=0.Ci 
RETURN 

4  PI=3.  14  159265 
HAL  PI =  P  I  *0 . 5 
TWP  I=2.0*Pl 
RTW'PI  =S CRT FIT  RP  1 1 
RHAPI  =SQRTF(  HAL  PI  1 
IF (A-4.  15,6,6 

C  A  LESS  THAN  FOUR 

5  AS  0  =  A  *A 
AFGUR=ASO*ASO 
AS  I  X=  AFCUR* ASL 
Z=- AS  C«  0. 25 
EXP-=£XPF  ( Z 1 

CALL  I Z  ER  0  (  -  Z  »  P.  ES  S  Z  ) 

CALL  IGNE  t -Z , BESSO 1 

EM1P=RHAPI«EXP»{BESSZ+AS0«( BESSZ+BESSO 1 *0 . 5 ) 

EM3P=RHAP I* EXP* 1 1  A SQ  +  3 . 0 1  * D ESS Z+ l AFOUR  +  4. O A  SC) *  1 B ESS Z  +  8E S SG ) *0 . 5 1 
EM  5P=  REAP  l  *  EX  P*  ( (  AFOUR+1  1  .  C*  A  S(J  + 1  5 . 0  )  *G  E  S  S  Z  +  (ASIX+  1  2 . 0*  A  FCLR  +2  3 . 0  * 
1  AS'QIMDCSSZ  +  BESSO  )*C.5) 

EM2P= ASG+2 . 0 
EMLP=AFOUR+S.C*ASO+S. 0 
EM6P= AS 1X+ io.0*AF0LR+72.C*ASQ+43.0 
S I GSG=EM2P- ( EM  IP**  2 ) 

SIGMA=SQRTF{$IGSQ 1  / 

EN1P2*EK1P«*2 
EM  1 P  3  =  E  M 1  P*  EM  VP 2 
EH 1 P4  =  EM  1 P  2*  *  2 
EM l P5  =  EM  1 P  4  *EM 1 P 
EMI P6=EM1 P3**7 

EM 3= EM 3 P - 3 . 0 * LM1P*  EM2P  +  2 . 0*EH 1 P3 

EMU  =  EM  i:  P-  li  .  0*  LM  1 P*  EM  5  P  +  6  .  C  *E  Ml  P2*  EM2P  -  3 . 0  *  E  M  1  Fli 

fcM5=EM  5P- 5 . 0*  EM  IP  *  EM  4P4 1 0 . 0*  EM  I P2  *EM3P- 1 J . C*  LM 1 P3*EM2P  +  4 . 0  *  EM  1 F5 
EM 6= EM 6 P-6 . 0*  lM  IP  *EM5P+ 1  5 . 0*EM  1  P2*£K4  p-20 . 0  *EM  1  P3*  EM5  P + 1  5 . 0  *  EM  1  P4* 
1  EM  2  P-  5  •  0*  EM  1 P  6 
TERM1 =EM5/( S I  CM A*  S I GSQ 1 
TERK2=EM4/(SI CSC* *2) -3.0 

TERM5-EM5/ t  S I GSO**  2*  S IGMA )- 1 C . 0* EM3/ ( S I GM A* S I GSC ) 
TEKM4=EM6/{$IGSQ**3)- 15. C*EM 4/tSIGSQ** 21+30. G 


485C 
4  84  C 
43  50 
4860 
4  87  C 
4  88  0 
4890 
49CC 
49  1  C 
49  20 

49  30 
494C 
4950 
4960 
4970 

4  98  C 
4990 
5000 
5010 
5C2C 

50  30 

50  40 
5C5C 
5C6C 
5070 
5Q0C 
5C9C 
5100 
5UC 
5120 

5  130 
5140 
5150 

51  6C 
5  1  7  C 
5100 
5190 

52  GC 
5210 
5220 
5230 
5240 
5250 
52  6C 
527C 
5230 
5290 
5300 


u> 

O' 


GO  TG  7 

6  AS  Q=A  « A 
ACUliE  =  A*A  SC 
AF0UR=ASC**2 
AF  I VE= AFOUR*A 
AS  I  X= ACUBfc**  2 
AS  V EN  =  AS  I  X* A 

EMlP=A+0. S/A+l  .0/(  8.C*ACU6EJ+3.0/M6.0*AFIVE)+75.0/{128.0*ASVEM 
SI GMA- 1 .0-1 .0/(4.0*ASQ)-S.C/ (32. 0*  AFOLR  ) -  97 . C / £1 2  S . 0* AS  I  X  } 

T6RM1 = 1 .O/ACUnE+3 .C/AFIVL 
TER  U2-- 3 . 0/ AFOUR 
TERM3=C  .0 
TERM4=C.C 

7  WHY=( V-EM 1 P ) / SIGMA 
WHY2=WFY**2 
WHY3=KHY*WHY2 

WHY  4= WFY2  *«  2 

WHY 5= WHY « WHY 4 

CALL  FHI ( WHY*  W»  ERF  »CPK1 ) 

■  PA V=C  PF I- ( EXPF ( -WHY2*0 . 5 ) /RT  WP I ) * ( ( WHY2- 1.0) / 6 . 0 ®T ERM 1  +  ( WHY  3-3 . 0  * 
1  WHY  )/2»-  ,0*T£RK2*(  WHY  4-6.  C'*WHY2-»3.C  )  / 1  20 . 0  *  T  ER^3  +  ( WHY5  -  10 . 0  *WF  Y5 
2+15.0*WHY ) /72C . 0*  T  ERM4 ) 

RETURN 

ENf) 

SUBROUTINE  PHI  (WHY, W, ERF ,CPHI ) 

W= ARS  F ( WHY ) /SORTF( 2.0} 

CALL  ERROP.R  (  W  ,  FRF  ) 

I F ( WHY-0. 000 OLl ) 6  0  0 , 6C1, 6 GO 

600  CPHI -C .5* ( 1 .0+ ( WM Y / A 8 S F ( W HY ) *  E  R  F » ) 

RE1  URN 

601  CPHI=C.5 
END 

SUBROUTINE  ERRCRR  (  W ,  E  R  F  } 

P=  ( (  ( ( 0 . 00  322975* W-0 .000394  46  I  ft  W+Q . 0274354  9 ) *  ki  +  G. Q8S64027 }  ft  W  + 

1  0.141  1 2  82 1 ) * W+ 1 . 0 

P=P«P 

P  =  P*P  -  ’ 

P=pftP 

ERF= 1  .C-l  .0/P 

RETURN 

END 

SUBROUTINE  I2lRO( EX, PESSZ } 

BESS  Z  =1.0 
EXSG  *  EXft-EX 
CEM0;M=  1 .0 
XHUtf=  cASG/4  .  C- 
FMUK=1.C 
DO  13C  1=1,10 


53  1C 
5320 
5330 
5340 
535  C 
5  360 
5370 
5380 
5390 
5400 
5410 

54  20 
5450 
5  4  40 


5460 

5470 

54  SO 
5490 
5  5  00 
5510 
5520 

55  3  C 
5540 
5  550 
5560 
5570 
5530 
5  5  90 
56CC 

56  1  C 
5  620 
56  3  C 
5  640 
5650 
5660 
567C 

56  8  C 
5690 
5700 

57  i  C 
5720 
57  50 
6740 
5750 
5760 
57  70 
57  GC 


I 


n=  i 

FNUM=FNUM*XNUK 

ccnck=ceno”*d 

1  GO  FESSZ=BE$SZ+{ FKUM / ( DENOM* *2 1 ) 
END 

SUBROUTINE  IONE (EX, BESSD ) 

PESSO  =1.0 

EX  S0=  EX*E  X 

DENI =1.0 

DEN2=  !  .C 

XNUM=  EXSQ/4 . 0 

FMUM=1.C 

DO  101  1=1,10 

C  =  I 

FNUM  =  F  NU.V*  XNUM 
DEN  1=  LEN1 *D 
DEN2=CEN2* ( D  + I .0) 

101  BESSO=BESSO  4 TNUM / ( DENI «DEN2 1 
BESSO  =BESSC*  CX*0 . 5 
END 
END 


v 


5790 
5800 
5610 
5820 
5820 
5840 
586C 
5660 
53  70 
5880 
5850 
5900 
59  10 
552C 
5930 
5940 
5950 
5960 
5970 
5900 
5990 


